Possible broken inversion and time-reversal symmetry state of electrons in bilayer 

graphene 
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With the two-band continuum model, we study the broken inversion and time-reversal symmetry 
state of electrons with finite-range repulsive interactions in bilayer graphene. In the state, there are 
overlapped loop currents in each layer. With the analytical solution to the mean-field Ifamiltonian, 
we obtain the electronic spectra. The ground state is gapped. In the presence of the magnetic 
field B, the energy gap grows with increasing B, in excellent agreement with the experimental 
observation. Such an energy-gap behavior originates from the disappearance of a Landau level of n 
— and 1 states. The present result resolves explicitly the puzzle of the gap dependence of B. 

PACS numbers: 73.22.Pr,71.10.-w,71.27.-(-a 



Introduction. Bilayer graphene (BLG) has attracted 
much attention because of its potential application to 
new electronic devices^"— Since the flat energy bands are 
sensitive to the electron-electron interactions, the ground 
state of the interacting electrons in BLG is highly dif- 
ferent from the noninteracting picture. Theories have 
predicted various gapped broken symmetry states, such 
as a ferroelectric-layer asymmetric state^^^ or a quan- 
tum valley Hall (QVH) state^^ a layer-polarized antifer- 
romagnetic (AF) state^iiii^ a quantum anomalous Hall 
(QAH) state^i^iii a quantum spin Hall (QSH) state^^ 
a charge density wave state^^ a superconducting state 
in coexistence with antiferromagnetism (SAF)j^ and 
the gapless statesJ^ii^, Though the experimental obser- 
vations on the ground state are controversial (except 
one experiment has found a gapless statei^), several 



experiment! 



,20-23 



have provided the evidence for the exis- 



tence of the gapped state at the charge neutrality point. 
In particular, a recent experiment performed on high 
quality suspended BLG^ has found that the ground state 
is gapped and the gap grows with increasing magnetic 
field B as 



12^2 



(1) 



with Aq ~ 1 meV and a w 5.5 meVT^^. This gap be- 
havior cannot be explained by the existing microscopic 
theories because they show that the gap is weakly de- 
pendent of the magnetic field. The gap growing with the 
magnetic field is a puzzle. 

In this paper, based on Varma's loop-current idea de- 
veloped in studying the mechanism of superconductiv- 
ity in cuprates^ we study the broken inversion and 
time-reversal symmetry state (BITRSS) of electrons with 
finite-range repulsive interactions in BLG. With the 
mean-field approximation (MFA) to the interacting elec- 
trons in a two-band continuum model (2BCM), we ana- 
lytically solve the eigenstates of the electrons and prove 
Eq. ([T]) for the gap behavior in the presence of the mag- 
netic field. 

The 2BCM. We first briefiy review the 2BCM. We start 



with the lattice structure of a BLG as shown in Fig. 
1. The unit cell of BLG contains four atoms denoted 
as A, B on top layer, and A' and B' on bottom layer. 
The lattice constant defined as the distance between the 
nearest-neighbor (NN) A atoms is a w 2.4 A . The energy 
of intralayer NN [between A (A') and B (B')] and inter- 
layer NN (between B and A') electron hopping are t « 
2.9 eV and ti « 0.31 eV, respectively. The Hamiltonian 
describing the noninteracting electrons is given by 



(2) 



where C^^ 



(c 



cL<T'CA'fea) (in which 4 



Aka 



creates an electron at A sublattice with momentum k 
and spin u), and B^ is given by 



Bk = 



Ft 
F D 



with as the 2x2 zero matrix, and 
F = 

D = 



4 
ek 



ek 
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-ti 
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and with the momentum k = {kx,ky) (in unit of a =1) 
confined to the first Brillouin zone. Note that Ck vanishes 
at the Dirac points K — (47r/3,0) and K' — = 
(— 47r/3,0). At ±K, Ck ~ eo{±kx + iky) where eo = 
\/3t/2 and k = {k^, ky) is measured from the Dirac point 
K [K') and confined to a small region of K {K'). We 
hereafter use the unit of eo = 1 for energy. For the carrier 
concentration close to the charge neutrality point, we are 
concerned with the states close to the zero energy. Taking 
the transformation Cka ~ TkVk with 



Tk 



1 1 

-D-^F 1 
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FIG. 1: (color online) Top view of BLG. Atoms A (A') and 
B (B') are on the top (bottom) layer. The parameters t and 
ti axe the electron hopping energies between the NN atoms 
belonging to the same layer, and to the neighboring layer 
above or below, respectively. 

and 1 as the 2x2 identity matrix, we have 

n"^^" = ( o' I) . 
^ (A 't) ■ 

Clearly, hk describes the valence and conduction bands 
close to the zero energy, while the eigenvalue of D is ±ti . 
The states connected by D belong to the bands of overall 
energy separation ±ii from the zero energy. For low- 
carrier doping, hk is the only part to be considered, which 
is the 2BCM^^'^^ for noninteracting electrons. Within 
the 2BCM, the A and B' sublattices are considered for 
top and bottom layers, respectively. 

The BITRSS. We here consider the interelectronic in- 
teraction effect. The interaction Hamiltonian is given as 

iJ' = ^ ^ vu,Vj5nu5nvj (4) 

ijll' 

where Snu is the electron number deviation from the av- 
erage occupation at site i on layer I (hereafter denoted 
as li for short), and vu^ij is the interaction between elec- 
trons at sites li and I'j. As long as the exchange effect 
is considered, the exchange interactions between elec- 
trons are finite-ranged (or short-ranged) because the bare 
Coulomb interactions are screened due to the electronic 
charge density fluctuations.— The operator product in 
the interaction Hamiltonian can be written as 

where is the density operator of spin ct, and J^^ j,^- — 

{(^UaCiija — Ciij„ciia)/2i IS proportional to the current op- 
erator. Therefore, the electrons can be coupled through 
the currents. Using the MFA and considering only the 
current couplings, we get 

H' = 2tJ2 ^u.i'iiJh'Acvj.- (5) 

IV ij(T 



FIG. 2: (color online) Current (Jfi ij) — ±Ji between site i 
and its NNs j on a layer. The + (-) sign is for the electron 
motion along (opposite to) an arrow. 

Here, ( J;"^ j,^) is a statistical average to be self- 
consistently determined with the mean-field theory. For 
simply exploring the physics, in this work, we consider 
the case that the average is nonvanishing only for the 
intralayer bonds and denote the average as {JUj/j} = 

±J/(|j — i\)Su'. Here, Ji{\j — i\) means a bond current 
between the two sites li and Ij. The sign factor ± de- 
pends on the direction of the electron motion from site 
i to site j. Under the assumption of broken inversion 
symmetry, we have Ji (of top layer) = - J2 (of bottom 
layer) . 

For the sake of illustration, we first consider the sim- 
plest case that there exist only the NN bond currents. We 
denote the NN interaction as vh^ij = v and the NN bond 
current as J;(|j — i|) = J/. Figure 2 shows the picture 
of the bond currents between the site i and its NN sites. 
The sign -I- (-) is for the motion along (opposite to) an 
arrow. The total current density passing through the site 
i is zero. For every site on a layer, one can depict such a 
figure. On the triangle lattice of the top or bottom layer, 
one then can find the triangle current loops surrounding 
every three NN sites. Therefore, Fig. 2 shows an equiva- 
lent view of the loop currents on a layer. To proceed, we 
rewrite Eq. ([S]) for the present case in momentum space, 

i?' = ^Az(fc)cLQfe., (6) 

kla 

with Ai{k) — J/u[8sin(/cj;/2) cos(\/3fci;/2) — 4sinfc2;] = 
Jiv^{k) [£,{k) is so defined for later use]. Since the low- 
energy states are considered as aforementioned, we take 
the leading term of the expansion of A;(fc) at the Dirac 
points ±_fir, 

A,(±ii:) = ±6%/3J;t; = s„s/A (7) 

where = ±1 for valley iiiT, s; = 1 (-1) for top (bottom) 
layer, and A = 6-\/3^if. The quantity A is the order 
parameter of the BITRSS. The signs Sy and s; refiect 
the broking of time-reversal and inversion symmetries, 
respectively. Equation ([7]) is an important result. 
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Now the effective 2BCM Hamiltonian is obtained as 



-SyA 



(8) 



with e„fc 



{sykx + ikyY/ti. The eigenvahies are 



|2 ^ 



±Ek and the corresponding eigen- 



wave-functions V'ufc are given by 



±R. 



^vk'-- 



vk 



where R^^. = ^/l ± s^A/Ek/V^, and (j)vk = arg(e„fc). 
For the system at the charge neutraUty point, the order 
parameter A is determined by 



27V 



9v 



{i('Aka('Aka)K' ' (c^Aka^'Aka) 



K 



g^SwA 

2V 



T^{fk-fi)lEk 



(9) 




where the first hne is the definition with the k summa- 
tion runs over the first Brillouin zone [with N the to- 
tal number of A (B') atoms on top (bottom) layer], the 
summation in the second line is separated into two parts 
over valleys K and K' with ^(fc) « S,{±K) = ±6^/3 be- 
ing used, and in the last line the averages are carried 
out using the eigen-wave-functions given above. Here 
= f{±Ek) with f{E) as the Fermi distribution func- 
tion, and V = \/3N/2 is the single-layer area (in unit of 
a — 1). For the ground state, we can integrate out Eq. 
(HD and get 



(10) 



where kc ~ ti is the momentum cutoff. 

To reproduce the experimental data |A| = 1 meV from 
Eq. (O , w ~ 2eo is needed. Unfortunately, this value of v 
is too strong to acquire for electrons in graphene. Using 
a = 2.4 A and the dielectric constant of graphene ^ 4, 
the typical value of v is obtained as 0.6eo- The problem 
stems from considering only the NN bond currents. If we 
include the contribution from the currents between long 
bond sites, the strength of order parameter A = 1 meV 
can be obtained. To repeat the MFA for the case of finite- 
range interactions, we note that {J^ ij) = ±>/;(|j ~ *l) 
depends on the vector d = j — i- The sign factor ± of 
all the currents passing through the site i is illustrated in 
Fig. 3. In this case, there are various overlapped triangle 
current loops in each layer. By the MFA, we obtain the 
similar equation for determining the order parameter as 
Eq. ([7]) but with v replaced with 
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FIG. 3: (color online) Phase for the current (Jujj) = ±J;(jz — 
j\). The sign -I- (-) is for the electron motion from site i as 
the center point to a red (white) site j. 



where v{d) — vu^ij, and — 47r/3. Here the d sum- 
mation runs over the sites of A (or B') sublattice. If the 
summation runs over only the six NNs, VcS reduces to v. 
For the general case of finite-range repulsive interactions, 
one has Vcs > v. 

Landau states. In the presence of a magnetic field B 
applied perpendicularly to the sample plane, we take the 
Landau gauge for the vector potential, A = {0,Bx). 
In this gauge, the y component momentum ky is a 
good quantum number. Replacing the variable x and 
the operator kx = —i^x with the raising and lower- 
ing operators and a, ky + Ex = ■\/_B/2(a^ -I- a) and 
kx — i^J~Bj2(a^ — a), we get from Eq. dH) 



A —LOcOu 
-{jJrO^ —A 



t2 



(12) 



for the quasiparticles at valley K and 

1 —A — Wcfl^ 



(13) 



for the quasiparticles at valley K' , and Wc — 2B/ti is the 
cyclotron frequency. The eigenvalue e„ and eigenvector 
ipn of hK, for n > 2, are given by 



Eji , 



Rn^n-2 
TP I I Rn 



where E^ = ^ A^ + ujln{n - \) , Rt = ^1 ± A/E^/^/2, 
and 4>n is the nth level wave function of a harmonic os- 
cillator of mass TO = ti/2 and frequency centered at 
Xc = —ky/B. For n = and 1, we have 



9 ^ 

d 



v{d) s,iv?{Kxdx) 



(11) 



eo,i = A, Vo,! = 



90,1 
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While at valley if', the solutions are 

e„ ^ -En, ^„ = (^i^;-) , n>2 
eo,i = A, = . 

We note that there is no state of energy — A for n = 
and 1. The particle-hole symmetry is no longer valid in 
the presence of the magnetic field, in qualitative agree- 
ment with the experimental observation.— The energy 
gap is i;gap = eo - {-E2) if A > or £^2 - eo if A < 0. 
The gap can be expressed as 

£;gap = |A| + + A2 (14) 

which is exactly the same form as Eq. ([Ij. The factor 
a in Eq. ^ corresponds to a 2^/2 /ti (in units of 
£0 = 1 and lattice constant a = 1), which is independent 
of the strength of the interactions. From the absolute 
values of the parameters t, ti and a given above, we get 
a « 5.2 meVT~^, which is close to the experimental data 
5.5 meVT-i. 

Recently, Zhu et al. have studied the same problem 
by numerically diagonalizing a mean-field Hamiltonian 
on a lattice of finite size because their model cannot be 
analytically solved;^ They have also reached the same 
conclusion that the gap grows with increasing B. 

Note that the present effective two-band Hamiltonian 
does not apply to the edge of the system because where 
the coordination number of each site is different from that 
in the bulk and some current loops disappear. Therefore, 
the form of the hopping term in momentum space and 
the self-energy are all different from that in the bulk. 



The ground state observed by the experiment^ is in- 
sulating at the charge neutrality point without external 
electric and magnetic fields. Also, the gap can be closed 
by an electric field of either sign perpendicular to the 
graphene plane. The experimental observation thus rules 
out the SAFifi as well as the gapless states. ^'^-i? The QAH 
and QSH states are excluded because the conductance 
~ Ae^/h yielded by their edge states is not observed. 
The ferroelectric-layer asymmetric stated""— or the QVH 
state^" not only cannot produce the gap behavior with B, 
but also contradicts to the observation of no net charge 
polarization between the layers. The order parameter of 
the AF stat&iiii^ is weakly dependent of the magnetic 
field. Therefore, among the existing proposed states, 
the overlapped loop-current state is the most possible 
candidate for the ground state of electrons in BLG. 

Summary. We have studied the BITRSS of the elec- 
trons with finite-range repulsive interactions in BLG. 
With the 2BCM, we have analytically solved the eigen- 
states of the mean-field Hamiltonian and obtained the 
gapped ground state at the charge neutrality point. The 
order parameter is odd for both interchanges of valleys 
and layers. In the presence of the magnetic field applied 
perpendicularly to the layers, we have obtained the spec- 
tra of the quasiparticles. Because of the breaking of time- 
reversal symmetry, a Landau level oi n = and 1 states 
disappears, resulting in the large energy gap growing with 
increasing the magnetic field. The present result is in ex- 
cellent agreement with the experimental observationi^ 
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Grant No. 10834011, and the Robert A. Welch Foun- 
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